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1. INTRODUCTION 

Underwater communications are a key technology used in the ocean [1]-[3] for millions of years, in 
large oceans, some marine mammals such as dolphins and whales have used acoustic waves to communicate 
among themselves [4]. As of late, requests for underwater communications have increased because of the 
continued expansion of human actions in the underwater location. For example, observing the environment, 
toxic waste monitoring, exploring underwater, summation of scientific data, marine archaeology, offshore oil 
exploration, strategic observation [5]-[7], military issues, ocean exploration [8], and so on. 

Common information medium for underwater communication includes radio waves, optical waves, 
magnetic waves and acoustic waves [9], [10]. Among them, the most popular medium for long distance 
communication is the acoustic wave because of its unique features in comparison to other applicants [11]. 
Radio waves are not preferable at low frequencies, as these waves are not able to travel long distances [12]. 
As a result of much higher attenuation in the water, and is rarely used for underwater communication [13]. 
Optical waves are not so attenuated, but are affected by absorption, scattering [14] and also high ambient 
light levels limiting the transmission range [15]. So, acoustic waves become the best option for short range 
and long range underwater communication [12]. 
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Underwater noise is caused by three main factors: ambient or ocean background noise, vehicle noise 
and intermittent noise, including biological noise such as shrimps, ice cracking and rain [16]. Ambient noise 
may also be defined as a combination of various sources that are unable to uniquely identify [17]. Ambient 
noise seriously impairs acoustic communication [18]. Whereas ambient noise is frequently approached like 
Gaussian, in practice it is colored noise with a decaying power spectral density (PSD). This decay rate is 
about 18 (dB/decade). In contrast to ambient noise, intermittent noise frequently contains important Non- 
Gaussian elements [19]. Since the sound attenuation in the ocean depends upon the operating frequency, the 
ocean acts as a low pass filter (LPF) for background noise level. Therefore, PSD background noise is 
generally (1/(f)^n) when the noise has lesser power at higher frequencies and greater power at lower 
frequencies [20]. There are a number of articles showing which noise isn't following a Gaussian distribution, 
especially in underwater communication systems. Underwater acoustic noise (UWAN) is characterized by 
the actual PDF with an extended tail shape, an emphasis on impulse behavior because of the high incidence 
of largeness amplitude noise events [16], [21]. 


2. UNDERWATER ACOUSTIC NOISE 

There are four components to the UWAN: turbulence, shipping, wind and thermal noises. Each 
component has different frequency effects [22]. PSD is calculated as proportional to the operating frequency 
f (kHz) in (dB re micro Pa per Hz) [16], [21]. Turbulence noise Nt (f) is generated in a system because the 
internal underwater movement, such as breakage waves [23]. The equation used to calculate turbulence noise 
is [24], [25]. 


Nt(f) = —30log f +17 (1) 
Shipping noise is caused by heavy shipping in water [22]. It can be computed as [24], [26]. 
Ns (£5)= 20(s — 0.5) + 26 log f - 60 log (f + 3*10—2) + 40 (2) 


Where: s represents the shipping activity factor between O (light) and 1 (heavy). The formula used to 
calculate wind noise is [27]. 


Nv (£ w) = 7.5.Vw + 20 log f - 40 log (f + 4*10-—1) + 50 (3) 


Where: 
w represents the wind speed in (m/sec) and is the wind factor with a value between 0 and 10. Thermal noise 
is a noise produced by randomly moving water molecules. It can be computed as [23]. 


Nth(f) = 20log f —15 (4) 


The total power spectral density is the sum of thermal, wind, ship and turbulence and may be calculated using 
(5) [24], [28]. 


Nrotai A) = Sx SN) +Ns (E s) + Nw (E w) + Nen(f) (5) 


Figure 1 presents the experimental deep water, background (ambient) noise PSDs for various values 
of shipping activities (s = 0,0.5 and 1) and wind speed with about (4.16 m/sec, 8 knots). Each noise source is 
found to be predominant within particular frequency range. Nt(f ) is dominant at very low frequencies 
generally (0.1-10 Hz), while Ns (f, s ) is the main factor contributing to noise in the frequency range (10-200 
Hz). Nw (f, w) is dominated in the frequency range (200-100,000 Hz). For frequencies above (100 kHz) 
thermal noise is dominant [22]. 

Figure 2 shows the experimental deep water noise PSD for various wind speeds (w = 0, 5 and 10) 
m/sec and shipping activities (s = 0.5). At fixed frequency, it is observed that the noise power spectrum 
density increases with increasing wind speed. Because, as the frequency increases, the effects of wind noise 
decrease. As a result, the dominance of that noise over a system becomes negligible [23]. However, this 
region is used by most acoustic systems [29]. When thermal noise is large the data rate is very bad. Basically 
the frequency choice criterion depends on the thermal noise value, since thermal noise is entirely dependent 
on the frequency of the signal. And since the frequency increases the thermal noise also increases [23]. 

These noise sources are dependent upon weather and other environmental factors. According to 
Figure | and Figure 2. The noise power spectrum lies in the region between both curves (1/f) and (1/f^3). 
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As a result, PSD background noise is generally (1/(f)^n) when the noise has lesser power at higher 
frequencies and greater power at lower frequencies [20]. 

The UWAN in shallow water has more variation in time and location of deep waters. Consequently, 
modeling or predicting is more difficult. Shallow water environments include depths of up to 200 m, 
while deep water environments throughout the oceans and seas are considered to be at depths greater than 
2000 m [30]. 
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Figure 1. PSD noise in (dB re micro Pa per Hz) with W = (4.16 m/sec, 8 knots) and various shipping 
activities (s = 0,0.5 and 1) 
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Figure 2. PSD noise in (dB re micro Pa/Hz) with different wind speeds (w = 0, 5 and 10) m/sec and shipping 
activities (s = 0.5) 


3. STATISTICAL PROPERTIES AND ESTIMATION 

In contrast to deterministic signals, random signals have the best statistical properties. So, statistical 
properties are functions of autocorrelation, PSD and probability density function (PDF). Are explained in this 
document in details with equations. 


3.1. Auto-correlation function 

A random signal can be distinguished by the way in which it is valued at a given time T, is 
dependent upon a different time’s value (t, T). It is the autocorrelation function (ACF), which is the cross 
correlation between a signal and itself [30]. ACF can be found in (6). 


Ra (t, T) = E[x(T) x (t + T)] (6) 
Where: - 
E [] represents an expected operator. When the random process is wide sense stationary (WSS), ACF may be 
represented as (7). 


Rxx (t, T) = R(t) = E[x(t). x(t + T)] (7) 
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The function is dependent on relative time instead of just absolute T. Within the normalized 
situation, (0 < R,,(t) < 1) which zero indicates that there is no dependence and another indicates high 
dependence [31]. 


3.2. Power spectral density and estimation 

The power spectral density represents the frequency of random processes. PSD is estimated through 
estimating the power spectrum. Two types of power spectrum estimation (PSE) methods typically exist the 
first is parametric and the other is non-parametric. Contrary to parametric PSE, non-parametric PSE does not 
make any assumptions about the data generation process. The relationship of the ACF to the PSD is 
determined Based on the Wiener-Khinchine theorem [32]. 


Sal = Joi, Rax (e7 dt, jev—T 8 
Equivalently, 
Ryx(t) = ee Sex (Pelit df (9) 


Letting f = 0 in (8) We get: 

Sx(0) = J RxCt) dt (10) 
In the same way, letting t = 0 in (9) Produces. 

Ryx(0) = Joy, Sux(f) af (11) 


In (10) indicates that ACF is integrable if (S,, (0) is less than oo) as well as (11) indicates that a PSD Sxx(f) 
is integrable when (Rxx (0) is less than œ). Two equations are typical cases of conventionally color noise. 
Yet, this is various from color noise types AD that have the property (Sxx (0) equal œ) [31]. 


Five commonly used nonparametric spectrum estimation exists includes: the periodogram, the 
modified periodogram, Bartlett spectrum estimation, Welch spectrum estimation, and Blackman-Tukey 
spectrum estimation. The estimated power spectrum of the periodogram is limited because the periodogram 
is not a consistent power spectrum estimate in which the estimated white noise variance does not decrease 
with the signal length. Other ways are trying to minimize that problem. The Welch method is usually 
considered the most consistent PSE method. Welch method divides data from time series into segments 
(which may overlap) calculate a modified periodogram for each segment and calculate the average power 
spectrum density estimate. The average of the modified periograms tends to reduce the variance in the 
estimation of a single periodogram for all data recording. Segment overlap increases the number of segments, 
thereby reducing the variance. The signal is sampled by using a normalized sampling frequency, so Welch 
PSE is [32]: 


k-1 
Tee ip) |? 
Preter(e W) = oe |w(n)X(n + iD)=/"™| (12) 
where: 

K: NO. of segments, L: the length of each segment, D: the offset between two successive segments, L-D: 
NO. of overlapping points and U is the normalization factor. The standardization factor may be expressed as 
(13): 


1 L-1 
U=} yw)? (13) 


3.3. Probability density function 

The standard noise model is Gaussian, additive, pixel independent and signal strength independent 
[33]. For many applications, Gaussian noise is mostly used as additive white noise for generating additive 
white Gaussian noise. This led to the probability distribution function given by: 


1 = 
P(z) = aE e 202 (14) 
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where -u: the mean value and o: the standard deviation. 

According to a number of publications, the UWAN does not follow the Gaussian (normal) 
distribution. Rather, UWAN is characterized by the with an extended tail shape reflecting an emphasis on 
impulse behavior because of the high incidence of largeness amplitude noise events [34]-[36]. UWAN is in 
accordance with the alpha-stable distribution class in which the characteristic equation is reversed. If we do 
not know the closed PDF format, the only solution is to use digital methods. Another modelling method 
involves empirically analyzing UWAN samples are taken from the seas and oceans [35] 

Student t distribution associated with Gaussian distribution which is characterized with wide tails. 
Student t distribution is applied to calculate the population mean for a few samples and the standard deviation 
of the population is unknown. The student’s t pdf is given by (15). 


v &-y)? 
E es 
p(x) m mel v ] 2 (15) 


Where: IT (-) represents the gamma function, u represents the location parameter, o represents the 
scaling parameter where (o > 0), and v represents the shape parameter where (v > 0), as v 
increases, the t-distribution of the student is approaching a Gaussian distribution as shown in 
Figure 3. 


Figure 3. Probability gravity function of the t-distribution of student 


4. RESULTS AND DISCUSSIONS 

A number of noise samples were taken and collected from Lake Hamrin during experiences 
conducted at the, Diyala, Iraq. (Latitude: 34.07° N; longitude: 44. 5829° E) in 29 November 2021. 
The signals were received using a hydrophone with a range (20~20,000) Hz, model (Dolphin-EAR DE200 
series), at various depths start with 1 meter, 3 meter and 5 meter to the bottom of the lake where the depth 
was 9 meter. Temperature (T) about 20 °C and W was about (4.16 m/sec, 8 knots). This is illustrated by 
Figures 4 to 6. 


Figure 4. Trial location of the experiment Figure 5. Field experiments in the Abu Dali district 
Tigris Beaches-Baghdad-Iraq, on 31 October 2020 
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Figure 6. Hydrophone model DolphinEAR DE200 series 
UWAN is registered by the hydrophone model (Dolphin-EAR DE200 series), that they convert into 
a discrete time representation for more treating and store in a computer at type (PC). We can use a modified 
periodic diagram technique from Welch to estimate the energy spectrum as the data are analyzed for different 
depths of 1 m, 3 m5 m and 7 m. Table 1 presents the details of the power spectrum estimation setup. 


Table 1. Factors specifications 


Factor Value 
Sample Frequency 8 KHz 
Type of window Hanning 
Number of points 2048 
Size of Window 256 
Segment Length 64 
Overlap percent 50% 


Figure 7 Displays the waveform time representation and autocorrelation function of data collected 
from the two different depths of Figure 7(a) 3 m shows time representation of the underwater noise at depth 
3meters, Figure 7(b) 5 m shows time representation of the underwater noise at depth Smeters, and Figure 7(c) 
7 m shows time representation of the underwater noise at depth 7meters and its clearly its more impulsive 
behavour than other depts. The number of noise samples is collected from Tigris River are considered 
correlated and colored noise due to the autocorrelation function of data collected is not like to the unit pulse 
This is different from white Gaussian noise. At t = 0, it represents the maximum magnitude of the noise 
power at ACF. Figure 8 describes the power spectrum estimation analyzed at various depths. It is obvious 
that UWAN has a decaying PSD. Whereas background (ambient) noise is frequently approached like 
Gaussian, in practice it is colored noise with a decaying PSD [19]. PSD background noise is generally 
(1/Cf )^n ) when the noise has lesser power at higher frequencies and greater power at lower frequencies 
[20]. However, the power spectrum density isn't uniform over the entire range of frequencies of interest as 
shown in the Figure 8, So the noise is colored, but not white. 

Figure 9. Demonstrates that the PSD of UWAN on a logarithmic scale with a deepness of 5 m is 
found from real field, where the rate of decay for the PSD of UWAN is about 20 (dB/decade). The actual 
field decay results compared to the [19] result differ by 2 (dB/decade), where the rate of decay for the PSD of 
ambient noise is about 18 (dB/decade). Figure 10. Indicates the PDF of UWAN found from the waveform 
time representation at a depth in Figure 10(a) 3 m shows pdf at depth 3m with degree of freedom 2.31, Figure 
10(b) 5 m shows pdf at depth 5m with degree of freedom 2.82 and Figure 10(c) 7m shows pdf at depth 7m 
with degree of freedom 2.93 and at this depth its clearly the pdf of underwater noise fitting to non-Gaussian 
noise. Better t location-scale proper for PDF of sample noise with use the distribution fitting tool in the 
Matlab. The results of the comparison between the Gaussian distribution and the Student t distribution make 
it clear that the amplitude of the UWAN usually follows the t distribution of Student. In Table 2 indicates the 
degree of freedom for different depths 3, 5 and 7 meter. According to the Table 2, the average degree of 
freedom is approximately 2.687. It is observed that the values freedom Degrees are small which mean that 
UWAN does not follow the Gaussian (normal) distribution. Rather, UWAN is characterized by the with an 
extended tail shape reflecting an emphasis on impulse behavior because of the high incidence of largeness 
amplitude noise events [34]-[36]. The analysis demonstrated that the characteristics of the underwater 
acoustic noise do not similar the added white Gaussian noise. The probability density function of UWAN 
usually follows the t distribution of Student, unlike the supposition of a Gaussian probability density function 
suggested in a previous study [37]. 
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Figure 7. The waveform time representation and ACF of the UWAN for two various depths at (a) 3 m, 
(b) 5 m, and (c) 7m 
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Figure 8. Welch PSD noise estimation at a depth 1 m, 3 m, and 5 m 
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Figure 9. Welch PSD noise estimation at a depth of 5 m logarithmically 
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Figure 10. Comparison between the student’s t -distribution and the Gaussian distribution for the probability 
density function of the real field samples noise at depths (a) 3 meter, (b) 5 meter, and (c) 7 meter 


Table 2. Freedom degrees for various depths 


Different depths (m) Freedom degrees (v) 
3 2.31 
5 2.82 
7 2.93 
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5. CONCLUSION 

The characteristics of the UWAN in Tigris River are particularly impulsive. PSD of the sample 
noise is taken from Lake Hamrin, Diyala, Iraq is not uniform over the entire frequency range and ACF of the 
sample noise isn't a delta function. Therefore, UWAN in Tigris River is considered colorful noise and Non- 
Gaussian distribution. Field trials presented that the noise power decreases from the surface of the lake to a 
depth of 5 meter, with a decay rate of about (20 dB/decade). The amplitude of the UWAN usually follows the 
t distribution of Student that is associated with the Gaussian distribution, but the tail's bigger. Consequently, 
UWAN does not take into account the assumption of white noise, meaning that Gaussian distribution is not 
applicable in the UWAN. 
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